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0^ ■ ABSTRACT 

^: 

^ ' Resonant relaxation is a novel form of two-body relaxation that arises in 

nearly Keplerian disks such as protoplanetary disks. Resonant relaxation does 
^ I not affect the semimajor axes of the particles, but enhances relaxation of particle 

eccentricities and inclinations. The equilibrium state after resonant relaxation 



> 



is a Rayleigh distribution, with the mean-square eccentricity and inclination in- 



'sj- I versely proportional to mass. The rate of resonant relaxation depends strongly 
on the precession rate of the disk. If the precession due to the disk's self-gravity 

I is small compared to the total precession, then the relaxation is concentrated 

(yQ \ near the secular resonance between each pair of interacting bodies; on the other 

^ I hand if the precession rate is dominated by the disk's self-gravity then relaxation 

' occurs through coupling to the large-scale low-frequency m = 1 normal modes 

^ I of the disk. Depending on the disk properties, resonant relaxation may be either 

I stronger or weaker than the usual non-resonant relaxation. 



1. Introduction 

The formation of planets from a disk of planetesimals is largely determined by two closely 
related processes: physical collisions and gravitational relaxation (e.g. ^tewart and Wetherill] 



1988 , Lissauer and Stewart 1993 ). Collisions result from close two-body encounters and drive 
the evolution of the mass spectrum of the planetesimals, while relaxation arises mainly from 
distant two-body encounters and drives the evolution of their phase-space distribution. The 
relative rates of the two processes are determined by the Safronov number, 6 ~ (fe/fr)^, 
where Ve is the surface escape speed and Vr is the rms random velocity of the planetesimals; 
relaxation dominates when ^ 1. Gravitational relaxation has two main effects on the 
random velocity of a particle in a disk: scattering drives a stochastic random walk towards 
larger values of Vr, while dynamical friction damps Vr- In general, scattering dominates 
for small particles while dynamical friction dominates for massive ones, as expected by 
equipartition. 
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The usual derivation of the rate of gravitational relaxation treats the evolution of the 
particle orbit as a sequence of uncorrelated two-body encounters with other particles. Such 
derivations neglect the effect of the gravitational field of the central star, but should be 
accurate in the early stages of planetesimal accumulation, so long as Vr ^ nrui where n is 
the mean motion and th = r{m/ M^^Y^^ is the Hill radius of the particle; here m and are 
the masses of the particle and the central star, and r is the orbital radius (|Ida 19901) . 



This papers investigates a qualitatively different type of gravitational relaxation, reso- 
nant relaxation, which we have already investigated in the context of spherical stellar systems 
( [Rauch and 'Iremaine 1996| ). Resonant relaxation arises in addition to the usual ("non- 
resonant") relaxation discussed above, when the potential in which the particles orbit is 
nearly Keplerian. To introduce the concept, let us imagine a time-exposure of a particulate 
Keplerian disk, in which the exposure time is much longer than the orbital period torb but 
shorter than the characteristic precession time tprec- Each particle appears in the photo- 
graph as an eccentric ellipse or wire. These wires exert torques on one another which remain 
roughly constant for ~ tprec! however, after several times tprec the configuration of the wires 
and the resulting torques will be quite different. These torques cause the angular momen- 
tum vectors of the particles — and thus their eccentricities and inclinations — to random walk, 
with the duration of one "step" in the walk being roughly tprec- However, the energies or 
semimajor axes are unaffected because the potential from the wires is nearly stationary. A 
closely related process is resonant friction, which damps the eccentricities and inclinations of 
massive particles without affecting their semimajor axes. The combined effects of resonant 
relaxation and resonant friction drive the phase-space distribution of particles towards the 
maximum-entropy state consistent with fixed particle energies and total angular momen- 
tum. The main goal of the paper is to estimate the resonant rates of excitation and damping 
of protoplanet eccentricities, to compare these briefly with the non-resonant rates, and to 
demonstrate that in some cases resonant relaxation dominates the eccentricity evolution. 

Resonant relaxation is a close cousin of secular perturbation theory in celestial mechan- 
ics, which likewise averages the Hamiltonian over times much longer than torb- Thus an 
alternative name for the process we are examining would be the oxymoron "secular relax- 
ation" . 

Nearly Keplerian disks can be parametrized by two dimensionless numbers 

^/ N q(r)n(r)r ^, ^ cr^ , ^ 

S{r) = - ) ' ^ = i , Q{r) = 1.07- ■ 



n{r) nr'^Ti(r) 7rG'S(r) ' n(r)r TTr^Ti(r)^ 

here g = w is the apsidal precession rate, S is the surface density of the disk, and is the 
rms radial velocity, related to the mean-square eccentricity by cr^ = |r^n^(e^). The second 
number is Toomre's Q-parameter ( |Binney and 'Iremaine 1987| ) and Q > 1 is required for 
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axisymmetric stability. If the precession rate is determined by the self-gravity of the disk, 
then generally S" ~ 1. For example, a disk with a power-law surface density S(r) oc has 

r(2 - i/?)r(i + i/?) 
^"^ r(| - |/?)r(i/?) ' "</^<^' 

in this case S" is independent of radius and varies only between 1 and 1.0942 in the range 
1<P<2. 

More generally, if there are other sources of precession then 15*1 ^ 1. In the rings of 
Saturn and Uranus, where the precession is dominated by the planetary quadrupole moment, 
S ~ —10^. In protoplanetary disks containing a thin disk of planetesimals embedded in a 
thicker disk of gas, S will be roughly the ratio of gas mass to planetesimal mass, S ~ 10^. 
We shall find that resonant relaxation is quite different in disks with 5* ~ 1 and 15*1 ^ 1. 

To simplify the discussion, we shall focus on the effects of resonant relaxation on plan- 
ets, that is, on bodies much larger than the planetesimals with which they are interacting. 
Also, we shall concentrate on resonant relaxation of the eccentricities, and generally ignore 
inclination relaxation except for a brief discussion in 

The equilibrium distribution of the particle orbital elements in a resonantly relaxed disk 
is determined in Estimating the rate of resonant relaxation is a more challenging task. 
A variety of tools can be used for this purpose, depending on the properties of the disk. We 
shall employ secular perturbation theory in to provide analytic estimates of the rate of 
resonant relaxation and resonant friction in disks. We carry out similar calculations using 
WKB density-wave theory in and numerical normal-mode calculations in The results 
are summarized and discussed in S|6|. 



2. Thermodynamic equilibrium 

The equilibrium eccentricity and inclination distribution of the particles in a resonantly 
relaxed disk can be determined by the methods of statistical mechanics. On timescales 
longer than the resonant relaxation time, the disk should be in the maximum-entropy state 
consistent with its total angular momentum and energy and mass distribution. Thus the 
equilibrium phase-space distribution function is ( [Ranch and Iremaine 1996| ) 

/(r, V, m) = m) exp(mb ■ J), (3) 

where E = — GM^,/r, J = r x v, and b and w{E,m) are determined implicitly by the 
total angular momentum and the distribution of particles in mass and energy, both of which 
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are conserved^. We orient the coordinate system so that b = be^ and introduce the actions 
(Jc, J, J,), where = {GM^af^ is the specific angular momentum of a circular orbit with 
semimajor axis a, J = Jc(l — e^)^/^ is the specific angular momentum, = J cos/ is the 
z-component of angular momentum, and e and / are the eccentricity and inclination. The 
conjugate angles are the mean anomaly, argument of perihelion, and longitude of ascending 
node (e.g. Prouwer and Clemence 1961|) . Then we have 

/(r, V, m) = w[E{Jc),m] exp^mbJ^), (5) 

where E{Jc) = —^{GM^,Y/J^. The advantage of this form is that the phase-space volume 
element is simply {2n)^dJcdJdJz, so the density of particles in action space is given by 

dN = {2nYw[E{Jc), m] exp{mbJz)dJcdJdJz. (6) 

If the eccentricities and inclinations are small we may approximate J and Jz by Jc(l ~ 
and Jc(l — — i/^), so 

dN = W{Jc, m) exp[-im6L(e2 + P)]dJcde^dP, (7) 

where W{Jc,m) = 2'7i'^w[E{Jc),m]J^ exp{mbJc)- This is the Rayleigh distribution (e.g. [Lis- 



sauer and Stewart 1993| ), with the further constraint that the mean-square eccentricity and 
inclination at different semimajor axes and masses are related by (e^), (P) oc l/(ma^/^). 



3. Secular perturbation theory 

We examine the relaxation of a planet of mass mp, semimajor axis ap, mean motion np, 
eccentricity ep <^ 1 and longitude of periapsis zup. We shall also use the complex eccentricity 
zp = epexp(zcup). The disk is axisymmetric and composed of planetesimals with masses 
rrij and orbital elements aj, nj, Cj, Wj, Zj, j = 1,...,N. The planet is assumed to be 
much more massive than the planetesimals, rrij <^ mp <^ M^. Kepler's law states that 

We assume that the eccentricities are small in the sense that ep,ej <^ \aj — ap\/ap. 
Then the secular evolution of the complex eccentricity is described by Lagrange's equations 



^In contrast, the equilibrium distribution function for non-resonant relaxation is 

f{r,v,m)=w{ni)exp{—ni(3E + mh-3), (4) 

where (3 and b are constants. This cannot be achieved in a Keplerian potential, since it requires that the 
density diverges as exp{GMi,m(3/r) as r ^ 0. 
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N 

■kjZj. 



( |Brouwer and Clemence 1961|) : 

N 

zp = igpzp - i ApjZj, 

J V 

4 = ig{ak)zk-iAkpZp-iJ2A 
Here 

The function g{ak) = and gp = wp are the rates of apsidal precession of the disk and the 
planet; we defer a discussion of these until §3 .31 . The Laplace coefficient 6^™^ (a) is defined 



by 

(1 — 2a; cos (j) + 



b^Ha) = ^f \^ (10) 



and we will use the result that 



(") - n ^ ^2 asa-^1. (11) 



3.1. Resonant relaxation 

Relaxation arises from the stochastic forces exerted on the planet by the planetesimal orbits. 
We concentrate first on disks in which the surface density is low, in the sense that IS"! ^ 1 
(cf. eq. m). In this case we can drop the term —iJ^jjtkA-kjZj which represents mutual 
perturbations among planetesimals from the second of equations (||). We also drop the 
term —lAkpZp which represents the forces from the planet on the planetesimals. Thus 
the planetesimal eccentricities satisfy Zk{t) = Zk{0) exp[ig{ak)t], where Zk{0) is the initial 
eccentricity. Substituting this result into the first of equations (^, and assuming that the 
initial planet eccentricity zp{0) = 0, we obtain 

zpit) = Y.ApM0) "''^''^"'^'^7''^^''''^ - (12) 
j 9p - 9{aj) 

We now compute the rate of change of the mean-square planetary eccentricity (cp), 
where (•) denotes an ensemble average. For axisymmetric disks the longitude of periapsis is 
uniformly distributed, so (2;* (0)^^(0)) = if j 7^ A;. Using this result we find 

|(4) ^ U'^) H- ^ .i:.i,(i,,o)i^)fcMi. ,13) 
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At large times we may use the relation lim^^oo sin(xt)/x = ttS{x), where S{x) is the Dirac 
delta-function; we also replace the sum over particles by an integral over semimajor axis. 
Thus 

where a = min(a, ap)/ max(a, ap), S(a) is the surface density of the disk, and (m) and (m^) 
are averages over the mass distribution. In this approximation, relaxation is entirely due to 
particles located at the secular resonance a<j where g{as) = gp. 



3.2. Resonant friction 



To evaluate the resonant friction on the planet, we return to Lagrange's equations (§) and 
once again drop the term —iJ2j^k^kjZj representing collective effects within the protoplan- 
etary disk. Then if we assume that zp,Zj oc exp^iut) we obtain the dispersion relation 



9p 



UJ 



E 



ApjAjp 



0. 



g[aj 



UJ 



(15) 



We replace the sum by an integral to obtain 

TT mp f J](a)ada o 



4/2 (") 



ma) 



0. 



9W 



LU 



(16) 



The subscript "L" is a reminder that the integral for Im(co') > is the analytic continuation 
of the integral defined for Im(ti;) < (the Landau contour). Equation (|TB]) can be solved 
directly for the eigenfrequency of the coupled planet-planetesimal system, but it is easier 
and more informative to focus on the case where the coupling represented by the integral 
in equation (|16D is weak enough that Re(ti;) ^ gp. Then the rate of eccentricity growth or 
damping from resonant friction is given by 



dcp 
~~dt 



— 2ep Im {uj) 



TT 9 mp^ r Tj(a)ada 

-epUp Im / a'^ 

4 ^ ^M^ Jl 



ma) 



g[a) - UJ 



;i7) 



Since the coupling is weak, lm{uj) is small. For Im(ci;) small and negative, we may use the 
relation 



lim 



sgn(e)7r5(x) 



to replace lm[{g 



UJ 



e^O e2 _|_ 

lmuj/[{g — gpY + (Imu;)^] by —T^5{g — gp). Thus we obtain 



de^p 



^ 2 



mp r ll{a)ada 



n{a)a^ 



(19) 
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and by Landau's prescription this is the correct expression for small values of lm{uj) of either 
sign. Thus resonant friction always damps the planet's eccentricity]^ 



Equation ([19| ) can be rewritten as 

~df ~ "T^^'^^M^ \dg/dloga\^ 
where all quantities are evaluated at the secular resonance a^. 



(20) 



The planetary eccentricity is in equilibrium when d{e^)/dt given by (IIID plus del/dt 
given by ([T9|) is zero. This occurs when 



^1/2 mp{m) 



(e^), (21) 



where the planetesimal parameters (m), (m^), a, (e^) are evaluated at the secular resonance. 
Equation (pT]) is closely related to the equipartition theorem (§0). 

There is an additional complication for a planetesimal disk embedded in a gaseous disk, 
as the gas disk may also have a secular resonance with the planet. The gas resonance will 
have a different location since the surface density of the gas disk is larger and its precession 
rate is modified by pressure. Equation (0) describes the resonant friction from the gas disk 
if S(a) is replaced by the gas surface density and g{a) by the precession rate of the gas 
elements. 



3.3. Precession rates 

The location of the secular resonance and hence the resonant relaxation rate is determined by 
the rate of apsidal precession of the planet and the disk particles. Precession can arise from 
several distinct sources. For example, the quadrupole moment of the central star induces 
precession at a rate 

= 2U ' ^^^^ 
where i?* and J2 are the radius and dynamical oblateness of the star. The quadrupole 
precession is negligible for most protoplanetary disks, although it dominates the precession 
in other systems such as planetary rings. 



^This result holds even for spherical systems, so long as the velocity distribution is predominantly tan- 
gential ( Ranch and Tremaine 1996 ). 
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The planet induces precession in the planetesimals. In Lagrange's theory the resulting 
precession rate of an object with semimajor axis and mean motion is 

2Gmp ab^^L(a) 
9{ak) = ^^N,p, N,,=N,, = - (23) 

Strictly this is the free precession, i.e. the precession that would occur if the planet orbit were 
circular. There is also forced precession due to the planet's eccentricity, which corresponds 
to the term —iAkpZp in the second of equations (||). 

The planetesimals also induce precession in the planet and each other. By an obvious 
extension of equation (PBD, the free precession of an object at semimajor axis a/^ is given by 

9i-^) = E ^AT,,. (24) 

j=i rikak 

However, the correct interpretation of this result is subtle for a continuous disk. The sub- 
tleties are illustrated by the observation that equation (|^) always predicts prograde free 
precession {g{a) > 0), whereas smooth, continuous, axisymmetric disks usually induce ret- 
rograde free precession (cf. eqs. ^ and 0). Moreover in the continuum limit (A^ —* oo, 
rrij ~ A^~^), the precession rate predicted by equation (p^ diverges as A^~^, whereas in a 
smooth continuous disk the precession rate is finite. 



The discrepancy arises because equation (|24D assumes that particle orbits do not cross 
(more precisely, that ej,ek <^ ~ o-kl/o-k for all j,k). This assumption is correct for an 
object with small free eccentricity that orbits outside the disk (or even within a gap in the 
disk), but incorrect if the orbit is embedded in a continuous disk, no matter how small the 
free eccentricity. Equation (^) can be applied to a continuous disk only if the complex 
eccentricity z{a) is a smooth function of semimajor axis, so orbits do not cross (i.e. if there 
are forced eccentricities but the free eccentricity is zero). 

The simplest fix for this shortcoming is to "soften" the gravitational potential of the 
particles, by replacing the Laplace coefficient ( p!OD by 



^"^"vrio (l-2«cos0 + a2 + e2)^' 

where e -C 1. In this case the expressions (|^) and (^) for Pj^ and Nj^ are no longer 
valid, since they were derived using recursion relations among unsoftened Laplace coefficients 
( [Brouwer and Clemence 1961| ); they must be replaced by 
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The precession rate predicted by equation p^ ) can then be apphed to a continuous disk 
as long as the spacing between adjacent particles is small compared to the softening, \aj — 
aj_i| <^ eaj. 

We shall also use e more generally, to indicate the typical dimensionless thickness or 
epicycle size in the disk; for |1 — a| ^ e the Lagrange equations (|) are not accurate. 

Yet another possible source of precession is the gravity from the thick gas disk that 
generally envelops the planetesimal disk in the early stages of planet formation. 



N. 



jk 



4max(afc, aj) 



3.4. The resonant relaxation rate 



To estimate the resonant relaxation rate we first determine the location of the secular reso- 
nance. We write the particle precession rate as g{a) = gQ^a) + gi{a) , where go is the smoothly 
varying precession rate due to the disk mass in gas and planetesimals, quadrupole moment of 
the central star, etc., while gi{a) is the precession due to the planet (eq. |2^). The precession 
rate of the planet itself is gp = go{cLp) (unless the disk properties change sharply at the 
planet's location, for example if there is a gap). We then expand go{a) in a Taylor series 
around ap; using equation (pH]) we then have 



ap 



1 mp 



na 



2tc Mi, dgo{a)/da 



1/3 



ap 



(27) 



Of course, this result is only valid if ap ^ \as — ap\ ^ eap] in other words the planet mass 
must be large enough that the secular resonance is separated by at least the typical epicycle 
size or disk thickness (~ eap). These constraints can be expressed approximately as 



Se' « ^ « S, 
rriD 



(28) 



where is the disk mass. The rates of resonant relaxation and friction are then given by 
equations (|l^ and (III 



1 d{el) (27r)4/3 (m2) Sa^ 



{ej>) dt 
1 dcp 



Cp dt 



-n 



(27r)^/3 Sa2 



1 dgo 



n d log a 

1/3 



1/3 



n d log a 



frnp\ 



-1/3 



(29) 
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where all quantities are evaluated at a^. Note that smaller planets relax faster, because the 
resonant radius is closer. 

These results have several limitations: (i) We have neglected the term —iYl,j^k^kjZj in 
equations (H), which describes collective effects arising from the self-gravity of the disk, (ii) 
The delta functions in equations ([T^) and ( |T9l) imply that as t — oo all of the relaxation or 
drag comes from particles located exactly at the secular resonance. More generally, at large 
but finite times the relaxation or damping is due to the few particles concentrated within an 
ever-narrowing resonant band Aa ~ \daj dg^jt] eventually there may be too few particles 
and too little mass in the resonant band for the formula to be valid, (iii) A related concern 
is that the eccentricities of the particles near resonance will become so large that nonlinear 
terms neglected in equations (§) will limit the strength of the resonant interactions. 

These concerns can be addressed using the complementary tool of density-wave theory, 
which treats the planetesimal disk as a continuous, coUisionless, self-gravitating fluid. 

4. Density-wave theory 

The disk is taken to have surface density radial velocity dispersion crj.(r), and epicycle 
frequency Kir) = n{r) — g{r), where as usual n(r) and g{r) are the mean motion and apsidal 
precession rate. We disturb the disk with a surface- density perturbation proportional to 
Re {exp[«(/'^ k{r)dr + m(f) — ujt)]}, and seek a dispersion relation that relates the frequency 
uj to the radial and azimuthal wavenumbers k and m. In the WKB or tight-winding limit 
\kr\ ^ 1 the dispersion relation for a coUisionless thin disk with a Rayleigh distribution of 
eccentricities is (e.g. [Binney and 'I'remaine 1957| ) 

K\r) - [mn{r) - u]' - 2nGE{r)\k\F f ~ f^^'^ ^4^] = 0; (30) 

where 

F{s,x) = -il-s')e-^tT^^.^ (31) 

and lj[x) is a modified Bessel function. 

For brevity we shall focus on the predictions of density-wave theory for resonant friction; 
most of our results apply equally well to resonant relaxation. Thus we are interested in 
density waves excited by a planet on an eccentric orbit. The slowly varying non-axisymmetric 
component of the planet's gravitational potential has the form (e.g. poldreich and Tremaine 
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1980) 



GmpCp I . d 



$(r, 0, t) = ^"7^" ( 1 - ":f ) ^i/2(«) cos[0 - U7p(t)], (32) 
2max(ap,rj \ da/ ' 



which corresponds to m = 1 and frequency uj = gp = Wp. Since \g\, \gp\ -C n in this case, 
we may simphfy equation (|30|) by neglecting terms of order (7^: 

2n{T)[gp - g{T)] - 27rGS(r)|fc|F (^-1, = 0, (33) 



or 

9p-9{r) ■ 



ar{r) 



X 



0, (34) 



where x = k a^/n . In the hmit x <^ 1 this simphfies to 



7rG'S(r)|A;| gp \kr\ 

9p - 9{r) = or — - = 1 - — -; 35 

n[r) g[r) b[r) 

In the WKB approximation, the planet excites a density wave satisfying the dispersion 
relation (p4D. The wave is excited near the secular resonance, where gp = g{r) and A; ~ 0. 
The wave propagates in the region where g{r) < gp; in the usual case where g{r) < and 
\g{r) \ decreases outwards, the wave propagates inwards as a leading spiral {k < 0); this is the 
long-wavelength branch of the dispersion relation. The wave becomes more and more tightly 
wound (I A; I increases) as it propagates. Eventually the wave encounters a turning point at 
X = 0.5841 (the maximum of the quantity in square brackets in eq. |3^) and thereafter 
propagates outward, continuing to become more tightly wound; this is the short-wavelength 
branch. The wave-propagation zone may be written 

0<9P-,(r)< 1.397!|M= 0.416^. (36) 

where Q is defined in equation (|T]). As the wave returns to the secular resonance \k\ ^ oo 
and the wave is damped by wave-particle resonances. 

The gravitational force from the leading spiral wave adds angular momentum — but al- 
most no energy, since the wave pattern is nearly stationary in an inertial frame — to the planet 
orbit. Thus the planetary eccentricity is damped. The damping rate can be determined from 
formulae in Goldreich and Tremaine (1980), and turns out to be precisely the same as equa- 
tion (|T^; this is a special case of the general result that the gravitational torques exerted 
on satellites by disks with a given surface-density distribution are largely independent of the 
nature of the collective effects in the disk ( [Goldreich and I'remaine 19791 ). 
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This analysis addresses several of the concerns raised at the end of the previous section. 
In particular, in the WKB approximation the rate of resonant relaxation and friction is not 
affected by the self-gravity of the disk; thus equations (|T^ and (|TU|) for the damping rate 
remain valid when collective effects are included. Moreover the collective effects ensure that 
a larger part of the disk is gravitationally coupled to the planet, not just the particles at the 
secular resonance: for example, equation (|3^) implies that most of the disk inside the secular 
resonance lies in the wave zone if Q <^ n/g; this is in contrast to waves in non-Keplerian 
disks or with m 7^ 1, where the much stronger condition Q ~ 1 is needed for a large wave 
zone. 

In fact, the principal limitation of the WKB analysis is that in many disks the collective 
effects are so strong that the WKB approximation is invalid. The "width" Aa^ of the secular 
resonance is roughly given by Xjl-Aa^ = 2n, where k is given by the WKB dispersion 
relation (^). We find 

Aa, = 2n ( ^) = 2nV' ( '-^) . (37) 

\n\dg/dr\J^^ V d^^Sr J 

The WKB analysis is only valid if the resonance width is much less than the distance between 
the planet and the secular resonance, which in turn requires 



mp 
vriD 



> S-^'\ (38) 



The failure of the WKB approximation is even more complete in disks with ~ 1. 
Since <C nr in thin disks, the dispersion relation ( P^ can only be satisfied if the quantity 
in square brackets is -C 1. This in turn requires that x -C 1 or x ^ 1; the latter condi- 
tion corresponds to the short-wavelength branch of the dispersion relation, which is rapidly 
damped by wave-particle interactions (Landau damping). Thus we focus on the case x <^ 1, 
where the dispersion relation is given by (^), which is the same as for a cold disk (eq. |33| 
with (Tr ^ 0). U Qp/g and S are of order unity, then (|35|) demands that \kr\ ~ 1 so the 
WKB approximation is not valid, however small the surface density of the disk may be. This 
situation contrasts with the case of a non-Keplerian potential or an azimuthal wavenumber 
m 7^ 1, where the WKB approximation \kr\ 3> 1 becomes more and more accurate as the 
surface density becomes smaller. Since the WKB approximation is invalid, we can only make 
an order-of-magnitude analytic estimate of the resonant damping rate. To do this, we set 



'E{as)a1 ~ mp) in equation (pO|), where mo is the disk mass. Since the precession is domi- 
nated by the disk self-gravity, we set \dg/dloga\ ~ nmp)/M^. Setting the Laplace coefficient 
and as to unity and neglecting other factors of order unity we obtain 



1 dei, mp 

~ — n — 

el dt M, 



n-f. (39) 
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The eccentricity damping rate is independent of the disk mass, except that when the disk 
mass becomes smaller than the planet mass the precession is dominated by the planet so 
that equation (^Sj) is no longer valid. 

For a more accurate determination of the rate of eccentricity damping in disks with 
~ 1, we must turn to numerical calculations. 



5. Excitation of normal modes 

The preceding discussion suggests that the large-scale, low-frequency, m = 1 normal modes 
of thin disks with 5^1 are largely independent of the velocity dispersion (eq. |35D. Thus 
we can determine the response of the disk to a planet using Lagrange's theory (eq. H) as a 
discrete approximation to a continuous disk with complex eccentricity z{a). 

We first determine the normal modes. If the terms involving the planet are dropped 
from the second of equations (^, we have 

N 

Zk = ig{ak)zk - i AkjZj, (40) 



where g{ak) is given by equation (p^ . This has solution 



^fc = E^exp(za;„t), (41) 

n=l 

where -ipk = {^kn-kY^'^o-k, bn is arbitrary, and tUn and c" = (c", . . . , c^) are the n^^ eigenvalue 
and eigenvector of the matrix M defined by 

Mk, = g{a,)6kj - ^(1 - 6kj)Akj; (42) 

that is, (M — u;„I)c"' = 0. The matrix M is symmetric and real, so all of the eigenvalues 
Un are real, and the eigenvectors c" can be chosen to be real and orthonormal, that is. 

Now suppose that a planet is present, with eccentricity zp oc exp(zti;t). Since the 
eigenvectors are complete, the forced response of the disk particles may be written in the 
form 

^ ft r" 
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solving the second of equations (|]) we find 

bn = — - — ^kclAkP. (44) 



1 ^ 



- ^ k=l 



Using this resuh to ehminate Zj from the first of equations (§) we obtain the dispersion 
relation for uj: 



N 1 N 



1 I 

gp-oo-Y: E ^Ap^A,pc]4 = 0. (45) 



This can be contrasted to equation ([TsD, which was derived by neglecting the self-gravity of 
the disk; in this case = Skn, = g{o,n), and equation (^5]) reduces to (^). 



As in the discussion following equation (|1^) we concentrate on the case where the planet 
mass mp is small compared to the disk mass nip,. We may then proceed as in equations 

(0) -(|T|) to obtain 

--^ = -27r E ^(^n - 9P) E j^Ap,A,pc]cl (46) 

n=l j,fc=l Y^i 

The eccentricity damping arises from resonances between the planet's precession frequency 
and the disk eigenfrequencies. In this derivation the spectral lines of the disk and planet 
have zero width, but in practice a number of mechanisms are likely to broaden these lines: 

(1) evolution of the planet's mass or semimajor axis during the planet-formation process; (ii) 
evolution of the disk mass or surface density due to infall or angular- momentum transport; 

(iii) damping of the disk normal modes from wave-particle resonances or other mechanisms; 

(iv) the planet's spectral line is broadened by 0(lma;) which is ~ gp{mp/mpi) for disks with 
S* ~ 1. Another route to a similar conclusion is to recall that the typical precession time 
g^^ is likely to be 10^-10^ yr (for a disk mass of 10^^-10^'^ Mq at a radius of 5 AU), and the 
likely formation time for the giant planets is 10^-10^ yr, during which the disk parameters 
and planetary orbit evolve substantially. Thus we expect the precessional spectra of the disk 
and planet to be broadened by at least of order 1%. 

To provide a crude representation of this effect — and to help ensure that the damping 
rates we obtain are well-behaved — we shall broaden the delta function into a Lorentz function 
by replacing 6{gp — Un) by T^~^hgp/[{gp — ojnY + h^g'pl^ a replacement which is exact in the 
limit — > +0; usually we pick h = 0.1. The use of a Lorentzian rather than some other 
broadening function is arbitrary but plausible. 

Figure 1 shows examples of damping rates computed by equation (^). We have exam- 
ined two disk models: the first has sharp edges, with surface density 

E (r) = i ^o{ro/ry-^ < r < tq 

^ 1 otherwise; ^ ^ 
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Fig. 1. — The damping rate d\og{ep)/dt for the eccentricity of a planet orbiting outside a 
Keplerian disk, as a function of the planet's semimajor axis. The curves labeled "1" are for a 
sharp-edged disk with surface density given by equation P7| ) and the curves labeled "2" are 
for a smooth disk (eq. ^). The solid curves include the effects of the disk self-gravity (eq. 
^61) and the dashed curves do not (eq. [T9|). The units are G = = tq = 1; the damping 
rate is proportional to the planet mass mp and is plotted for mp = 1. 
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the second is smooth and is obtained by multiplying Si(r) by a window function in logr: 

S2(r) = Si(r)sin^[7rlog(r/ri)/log(ro/ri)]. (48) 

We choose rj = O.lro, h = 0.1, and the results are expressed in units in which G = = 
tq = 1. The eccentricity damping rate is proportional to the planet mass mp and is plotted 
for mp = 1; in this approximation the damping rate is independent of the disk surface 
density Eq. The planet is assumed to orbit outside tq, and the damping rate is plotted as 
a function of the planet's orbital radius. The solid curves are obtained from equation 
which includes the effects of the disk's self-gravity, while the dashed curves do not (eq. 



The damping rate for the sharp-edged disk exhibits a resonant peak near an orbital 
radius of 1.2. The peak arises from a resonance between the planet's precession rate and 
the fundamental m = 1 mode of the disk (i.e. the normal mode with zero nodes). At larger 
orbital radii the damping rate is determined by the tail of the Lorentz function associated 
with this and other normal modes (and thus depends on our arbitrary choice of a Lorentzian 
profile with width h = 0.1). 

The damping rate is lower for the smooth disk; in this case the eigenfrequencies of the 
disk are higher than the precession frequency of the planet, so the damping rate for a planet 
outside the disk is everywhere determined by the tails associated with each resonance. 



These results can be re-stated by writing the heuristic formula (|39|) as 

' -/^^; (49) 



1 de% mp 



-p 

then the dimensionless factor / is approximately the same as the ordinate of Figure 1 for 
planets orbiting at semimajor axes near the disk edge. The important feature of Figure 1 is 
not the numerical damping rates, which depend on our arbitrary broadening function, but 
that even modest broadening leads to substantial damping over a large radius range. 

Computing the eccentricity damping rate for a planet embedded in a disk is more 
complicated, since (i) the planetary precession rate and hence the damping will depend 
sensitively on whether there is an annular gap in the disk surrounding the planet; (ii) the 
planetary potential will have higher spatial frequencies and thus will interact with a richer 
set of disk normal modes; (iii) the approximation of small eccentricity is only valid if e, ep -C 
\a — ap\/ap, and this is much more difficult to satisfy for nearby particles; (iv) non-resonant 
friction due to close encounters between the planet and disk particles will also play an 
important role in eccentricity evolution (see §|| below). 
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6. Discussion 

We have investigated resonant gravitational relaxation, which arises in particle disks orbiting 
in a nearly Keplerian potential. Resonant relaxation does not affect semimajor axes, and 
tends to produce a Rayleigh distribution in eccentricity and inclination, with mean-square 
eccentricity or inclination inversely proportional to mass (eq. |^). 

The nature of resonant relaxation depends strongly on whether apsidal precession is 
dominated by the self-gravity of the disk, a condition which is parametrized by S{r) (eq. |1|). 
If other mechanisms dominate the precession (IS*! ^ 1) then resonant relaxation requires 
that a secular resonance is present in the disk, which in turn requires mp/niD <^ S, where 
nip and mp) are the planet and disk masses (eq. ^S]). If the resonance is present, the rate of 
damping of the planet's eccentricity from resonant friction is roughly (eq. ^) 

l^^^.n(!^Y's'r'('^y'\ (50) 

ep dt \MJ \M^J ' ^ ^ 

and the rate of stochastic excitation of the eccentricity is 

rip r-, 

(51) 



1 d{el) 




(-) 


1 dep 


(el) dt 




\mp/ 


ep dt 



These formulae also require that (i) the planet mass is large enough that the secular resonance 
is separated from the planet by more than the characteristic disk thickness or epicycle size 
(eq. |2^); (ii) the separation exceeds the characteristic width of the resonance caused by the 
disk self-gravity, mp/mp, ^ S~^^'^ (eq. P^). 

If on the other hand the apsidal precession is determined by the disk mass {S ~ 1), 
the planet interacts with the large-scale low-frequency normal modes of the disk rather than 
with particles near the secular resonance. Calculating the rate of resonant relaxation in this 
case is harder (cf. especially if the planet is embedded in the disk. For heuristic purposes 
the rate of eccentricity decay from resonant friction may be written as (eq. ^) 



1 de% nip 
~ — n — 

el dt M, 



^T7^- (52) 



A key point is that the damping rate is independent of the disk mass m^i so long as the disk 
mass determines the precession rate. 



A major uncertainty in equation ( ^21) is whether the planetary precession rate is in 
resonance with one of the disk eigenfrequencies; this usually requires either a continuous 
spectrum of eigenfrequencies or that dissipative or evolutionary processes broaden the dis- 
crete disk eigenfrequencies so that they overlap. Because the characteristic frequency of the 
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disk modes is of order n{mD/M^^) <C n, even slow evolutionary or damping processes cause 
considerable broadening. An important next step will be to understand the low-frequency, 
large-scale m = 1 eigenmodes of nearly Keplerian disks. 

For comparison, the non-resonant eccentricity damping rate for a planet embedded in a 
disk may be written as ( ptewart and Wetherill 1988| , P^issauer and Stewart 1993| ) 



1^^ mpm^/^y (53) 
el dt \vrJ 

where Vr is the rms random velocity of the planetesimals, Vc is the circular speed, and In A 
is the Coulomb logarithm. As described in this result is valid so long as Vr ^ nrn where 
vh is the planet's Hill radius ([Ida 1990|) . 

If a planet orbits at a distance Aa from a narrow ring, the non-resonant damping rate 
becomes ( [Goldreich and 'I'remaine 198U| , eq. 31) 



1 dcp nip nxD 




^^^^ ^^^^ 



el dt A't 

These crude formulae suggest that either resonant or non-resonant eccentricity damping 
can dominate in a given disk, depending on the disk mass, velocity dispersion, and other 
parameters. In typical protoplanetary disks, non-resonant damping probably dominates if 
the planet is embedded in the disk. However if there is a gap Aa/ap ^ {rrip,/ M^^y/^ then 
resonant damping may be stronger. 

Resonant inclination relaxation is similar to resonant eccentricity relaxation. The princi- 
pal difference is that the precession due to the self-gravity of the disk is g ^ nmp, I M^^ivcl ^ 
stronger by a factor of order Vc/v^- Thus the minimum value of the parameter IS"! ~ v^/vj. ^ 
1. 

Resonant relaxation is a complex process whose efficiency depends sensitively on proper- 
ties such as the precession rate, the presence of gaps around planets, etc. (although it is much 
less sensitive to the rms random velocity in the disk than non- resonant relaxation). Moreover 
many of the tools that are usually employed to investigate disk relaxation (two-body scatter- 
ing, Fokker-Planck approximation. Hill's problem, etc.) are too crude to describe resonant 
relaxation. Nevertheless this process is likely to play an important role in the evolution of 
a variety of nearly Keplerian disks, including both protoplanetary disks and other objects 
such as stellar disks surrounding black holes at the centers of galaxies. 
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